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We report on the status of our QCD thermodynamics project. It is performed on 
the QCDOC machine at Brookhaven National Laboratory and the APEnext machine at 
Bielefeld University. Using a 2+1 flavor formulation of QCD at almost realistic quark 
masses we calculated several thermodynamical quantities. In this proceeding we show the 
susceptibilites of the chiral condensate and the Polyakov loop, the static quark potential 
^ ! and the spatial string tension. 

O 

1. Introduction and Lattice Setup 

The calculation of QCD thermodynamics from first principle is important for various 
research areas such as Heavy Ion Phenomenology, Cosmology and Astrophysics. Lattice 
QCD enables us to carry out such calculations. Especially for HIC phenomenology it 
is mandatory to improve estimates on some basic thermodynamic quantities which have 
been obtained in previous lattice calculations. Since thermodynamics of lattice QCD re- 
quires huge computational resources, it is difficult to perform an ideal simulation. Recent 
studies tell us that quark masses and the number of flavors strongly affect thermodynamic 
quantities [[!]. Reliable continuum extrapolations are of tremendous importance as well 
[12]. Therefore, it is our goal to study QCD thermodynamics with almost realistic quark 
masses on the QCDOC machine at Brookhaven National Laboratory and the APEnext 
machine at Bielefeld University. The calculation is performed with Nf = 2 + 1, which 
means 2 degenerate light quarks and one heavier quark on lattices with N t = 4 and 6. 
The lightest quark masses of our simulation yields a pion mass of about 150 MeV and 
kaon mass of about 500 MeV. 

For such calculations we adopt the p4fat3 quark action, which is an improved Staggered 
quark action [ 3 , with a tree- level improved Symanzik gauge action. By using the p4fat3 
action, the free quark dispersion relation has the continuum form up to 0(p A ), and taste 
symmetry breaking is suppressed by a 3-link fattening term. The action also improves 
bulk thermodynamical quantities in the high temperature limit [ 3J. The improvements 
are essential to control the continuum extrapolation on rather coarse lattices, i.e. iVj = 4 
and 6. The gauge ensembles are generated by an exact RHMC algorithm [0]. 

As a status report of the project, in this proceeding, we present several thermodynamics 
quantities, which are susceptibilities of the light and strange quark chiral condensate, the 
Polyakov loop susceptibility, the static quark potential, and the spatial string tension. 
The details of the critical temperature calculation are given in our recent paper [|3]. 
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Figure 1. The susceptibility of the chiral condensate (left) and the susceptibility of the 
Polyakov loop (right) on 8 3 x 4 lattices. In the left panel filled and open symbols are results 
for light and strange quarks, respectively. The lines are calculated by the multi-histogram 
re-weighting technique. 



2. Order Parameters and Susceptibilities 

To investigate the QCD critical temperature and phase diagram, order parameter of 
the QCD transition are indispensable. In the chiral limit the chiral condensate (4>ip) 
is the order parameter for the spontaneous chiral symmetry breaking of QCD. On the 
other hand in the heavy quark limit the Polyakov loop (L) is the order parameter of the 
deconfinement phase transition. For finite quark masses, these observables remain good 
indicators for the (pseudo) critical point. Especially their susceptibilities are useful to 
determine the critical coupling (3 C in numerical simulations. 

Figure 1 shows the susceptibilities of the chiral condensate and the Polyakov loop. 
Their peak positions define the point of most drastic change of each order parameters, 
i.e. the (pseudo) critical point of the QCD transition. The results are interpolated in the 
coupling (3 by using the multi-histogram re- weighting technique 

The strength of the transition decreases with increasing quark masses, this is reflected in 
the decreasing peak height of the chiral susceptibilities. We calculate these susceptibilities 
on lattices with aspect ratios of N s /N t = 2 and 4. Since we see a rather small volume 
dependence the results suggest that the transition is in fact not a true phase transition 
in the thermodynamic sense but a rapid crossover. The peak position of the chiral and 
Polyakov loop susceptibilities are almost identical, i.e. the chiral and deconfinement 
transition occur at almost the same temperature. The discrepancy between the peak 
positions shrinks with increasing volume. 

3. Scale Setting and the Heavy Quark Potential 

The lattice scale is determined from the heavy quark potential V(r) which is extracted 
from Wilson loops. The Wilson loop expectation values are calculated on 16 3 x 32 lattices 
with APE smearing in spatial direction. The spatial path in a loop is determined by the 
Bresenham algorithm We calculate the string tension, a and Sommer scale ro, which 
is defined [ |7j as the distance where the corresponding force of the static quark potential 
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Figure 2. Heavy quark potentials scaled by ro for various cut-offs and quark masses 
(left) and the fit-range dependence of r (right). The latter shows the result for N t = 
4,m q /m s = 0.026/0.065 and fits with the 3-params. fit-form. 



matches a certain value suggested by phenomenology: r 2 ^ \ r=ro = 1.65. To remove short 
range lattice artifacts we use the improved distance, r imp , which is defined as 

d 3 k e ikr 



Anr imp J (27r) 3 4E J (sin 2 |) + isin 4 |' 

In our lattice setup, we find almost no mass and cutoff dependence in the potential 
scaled by r at N T = 4 and 6 (Fig.2(left)). As discussed in previous studies we also 
find no string breaking effects even at large r. 

To estimate systematic uncertainties of the potential fit, we performed several types of 
fits, e.g. different fit-ranges in r (see Fig.2(right)) and fit-forms (3 & 4 params. fits), 

V(r) =C + — + ar imp , V(r) = C+- + ar + d{— - -). (2) 

The differences in the mean values of the fits are evaluated as a systematic uncertainty 
of the scale setting. From a combined quark mass and cut off extrapolation of tqT c = 
(ra/a)(aT c ) we finally obtain a critical temperature of T c = 192(7) (4) MeV at the physical 
point [E]. Here we used tq = 0.469 fm to set the scale. The first error summarizes all 
statistical and systematic errors on r and the critical couplings j3 c and the second error 
reflects the remaining uncertainties in the extrapolation. 

4. Spatial String Tension 

Let us now discuss the calculation of the spatial string tension which is important to 
verify the theoretical concept of dimensional reduction at high temperatures. The spatial 
string tension is extracted from the spatial static quark "potential" (from spatial Wilson 
loops). We use the same analysis technique as for the usual (temporal) static quark 
potential. 

At high temperature, the spatial string tension cr s {T) is expected to behave like 



l<j s (T) = cg\T)T. (3) 
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Figure 3. Temperature dependence of the spatial string tension for Nf — [ ITU] and 
iVj = 2 + 1. Both dotted lines are fits with Eq. |H1 Note, that the scale on the horizontal 
axis is logarithmic. 



Here g 2 (T) is the temperature dependent coupling constant from the 2- loop RG equation, 

9 -(T) = 2Mn£ + £ln(21n£). (4) 

If dimensional reduction works, the parameter "c" should be equal to the 3-dimensional 
string tension and should be flavor independent. 

Our 2+1 flavor result yields c = 0.587(41) and A CT /T C = 0.114(27), obtained by a fit 
with Eq. El On the other hand, we plot in Fig |3] also the quenched result [[TU] which gives 
c = 0.566(13) and A a /T c = 0.104(9). We thus find that the parameter "c" is - within 
statistical errors - independent on the number of dynamical flavors and that dimensional 
reduction works well even for T = 2T C . This analysis can and will be refined in the future 
by taking into account higher order corrections to Eq. |3] [lllj. 
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